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Strong interactions are addressed in connection to extracting \ V cb \ and determining heavy quark parameters. A comprehensive approach 
allows a robust analysis not relying on a l/m c expansion; a percent defendable accuracy in \V c b\ becomes realistic. Some of the heavy 
quark parameters are already accurately known. We have at least one nontrivial precision check of the OPE at the nonperturbative level 
in inclusive decays. The alleged controversy between theory and BaBar data on {M\) is argued to be an artefact of oversimplifying the 
OPE for high cuts in lepton energy. Consequences of the proximity to the 'BPS' limit are addressed and their accuracy qualified. It is 
suggested that theory- wise B— >£> £v near zero recoil offers an accurate way for measuring | | . 



Total semileptonic decay rate T S \(B) is now one of the best 
measured quantities in B physics. This opens an accu- 
rate way to extract the CKM angle \ V c i,\ where systematic 
uncertainties can be thoroughly studied. Nonperturbative 
effects are theoretically controlled by the QCD theorem 
[ [l] which established absence of the leading Aqcd i 'mt 
power corrections to total decay rates. It is based on 
OPE and applies to all sufficiently inclusive decay prob- 
abilities, semileptonic as well as nonleptonic or radiatively 
induced. The theorem relates the inclusive B widths to 
(short-distance) quark masses and expectation values of lo- 
cal fe-quark operators in actual B mesons. To order 1/wt? 
they are kinetic and chromomagnetic averages jj^, fA hav- 
ing a straightforward quantum-mechanical meaning. To 
order 1 /m? two new operators, Darwin and LS appear for 
b — > civ with the averages and p 3 LS , also analogous to 
those well known from atomic theory. The general expan- 
sion parameter for inclusive decays is energy release, in 
b — > civ it constitutes m\,-m c ^ 3.5 GeV. For a recent 
review, see Ref. [0. 

Within the OPE approach calculating everything suffi- 
ciently inclusive has become a technical matter (although 
not necessary simple as illustrated by perturbative correc- 
tions). Yet care must be taken to treat properly certain sub- 
tleties; some of them will be discussed. 

Heavy quark masses are full-fledged QCD parameters en- 
tering various hadronic processes. The expectation val- 
ues like /ij likewise enjoy the status of observables. The 
masses and relevant nonperturbative parameters can be de- 
termined from the B decay distributions themselves 
Nowadays this strategy is being implemented in a number 
of experimental studies. 

Some progress has been recently made on the theoretical 
side as well. A resummation of the dominant perturbative 
corrections in the context of Wilsonian OPE was performed 
for the semileptonic width and showed safe convergence [ 
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|5]|6), eliminating one of the major theoretical uncertainties. 
Higher-order nonperturbative corrections were studied and 
a new class of nonperturbative effects identified associated 
with the nonvanishing expectation values of the four-quark 
operators with charm field, (B\bccb\B) [ 6|. The so-called 
"BPS" regime was proposed as an approximation based on 
the hierarchy fi^-fi^ <sz ft^ suggested by experiment, and 
its immediate consequences analyzed [ 7|. The new results 
extending them will be addressed in the end. 

The new generation of data provides accurate measure- 
ments of many inclusive characteristics in B decays. At 
the same time, the proper theoretical formalism gradually 
finds its way into their analyses. To summarize the current 
stage: 

• We have an accurate and reliable determination of some 
heavy quark parameters directly from experiment 

• Extracting \V c b\ from T S \{B) has good accuracy and solid 
grounds 

• We have at least one nontrivial precision check of the 
OPE at the nonperturbative level. 

The latter comes through comparison of the average lepton 
energy and the invariant hadronic mass, as explained in the 
following section. 

Present theory allows to aim a percent accuracy in \V c b\- 
Such a precision becomes possible owing to a number 
of theoretical refinements. The low-scale running masses 
rribifi), m c (jj), the expectation values f£(fi), Hq(i-i), ■■■ are 
completely defined and can be determined from experi- 
ment with an in principle unlimited accuracy. Violation 
of local duality potentially limiting theoretical predictabil- 
ity, has been scrutinized and found to be negligibly small 
in total semileptonic B widths [ 8 1; it can be controlled ex- 
perimentally with dedicated high-statistics measurements. 
Present-day perturbative technology makes computing a s - 
corrections to the Wilson coefficients of nonperturbative 
operators feasible. It is also understood how to treat higher- 
order power corrections in a way which renders them sup- 
pressed [|6|. 

High accuracy can be achieved in a comprehensive ap- 
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proach where many observables are measured in B decays 
to extract necessary 'theoretical' input parameters. It is 
crucial that here one can do without relying on charm mass 
expansion at all, i.e. do not assume charm quark to be re- 
ally heavy in strong interaction mass scale. For reliability 
of the 1 /m c expansion is questionable. Already in the 1 /Mq 
terms one has -4 > 14-W; even for the worst mass scale in 

m r mi 

c b 

the width expansion, -. — 1 —tt is at least 8 times smaller than 
\. There are indications [ 7 1 that the nonlocal correlators 

nic — 

affecting meson masses can be particularly large - a pattern 
independently observed in the 't Hooft model [ 9 1. This ex- 
pectation is supported by the pilot lattice study [ 10 1 which 
- if taken at face value - suggests a very large value of a 
particular combination pl K +pl entering in the conventional 
approach. On the other hand, non-local correlators are not 
measured in inclusive B decays. 

The approach which is totally free from relying on charm 
mass expansion [ II II was put forward at the previous 
CKM-2002 Workshop at CERN. It allows to utilize the full 
power of the comprehensive studies, and makes use of a 
few key facts [|4][3]: 

• Total width to order 1 /ml is affected by a single new 
Darwin operator; the moments also weakly depend on p 3 s . 

• No nonlocal correlators ever enter per se. 

• Deviations from the HQ limit in the expectation values 
are driven by the maximal mass scale, 2mj (and are addi- 
tionally suppressed by proximity to the BPS limit); they 
are negligible in practice. 

• Exact sum rules and inequalities which hold for properly 
defined Wilsonian parameters. 

The original motivation for the precision control of strong 
interaction effects in B decays was accurate determina- 
tion of \V c b\ and \V„b\ as means to cross check the Stan- 
dard Model. For theorists it is clear, however that interest- 
ing physics lies not only in the CKM matrix; knowledge 
of heavy quark masses and nonperturbative parameters of 
QCD is of high importance as well. 

Some of the HQ parameters like p c are known before- 
hand. Proper field-theoretic definition allows its accurate 
determination from the B*-B mass splitting: p c (l GeV) = 
0. 35^Q2 GeV 2 [ 7|. A priori less certain is p 2 n . How- 
ever, the inequality p 2 > p? G valid for any definition of ki- 
netic and chromomagnetic operators respecting the QCD 
commutation relation [D,-, D^] = —ig s Gjk, and the corre- 
sponding sum rules essentially limit its range: p 2 (\ GeV) = 
0.45 +0.1 GeV 2 . 

Lepton and hadron moments. Moments of the charged 
lepton energy in the semileptonic B decays are traditional 
observables to measure heavy quark parameters. Moments 
of the invariant hadronic mass squared M\ in semileptonic 
decays is another useful set of observables. Their utility 
follows from the observation [ 4| that, at least if charm 
were heavy enough the first, second and third moments 



would more or less directly yield A, p 2 and p 3 D . Precision 
measurements of the B — > X s + y decays also may yield a 
number of constraints. Details of the present situation can 
be found in the experimental review talks ni2lll3llT4lll5l . 

Let me briefly illustrate how this strategy works number- 
wise. Lepton energy moments, for instance, are given by 
the following approximate expressions (b — > u decays are 
neglected): 

(E c ) = 1.38 GeV+0.38[(m/,-4.6 GeV) -0.7(m c - 1.15 GeV)] 
+ 0.03(/4-0.4GeV 2 ) - 0.09(p^-0.12 GeV 3 ) , 

{(E c - (Ec)) 1 ) = 0.18 GeV 2 + 0. 1 [(m b -4.6 GeV) - 

0.6(m c -1.15 GeV)] + 0.045(yu 2 -0.4GeV 2 ) 
-OX^Oo^-O.^GeV 3 ) , 
{(Ec-{E{)f) = -0.033 GeV 3 - 0.03 [(m fe -4.6GeV) 

- 0.8(m c - 1 . 15 GeV)] + 0.024(yK 2 - 0.4 GeV 2 ) 

-0.035(^-0. 12 GeV 3 ) (1) 

(all dimensionful factors are given in the corresponding 
powers of GeV). The moments depend basically on one 
and the same combination of masses m& — 0.65m c ; depen- 
dence on p 2 is rather weak. To even larger extent this 
applies to the CLEO's cut moments R\, R2 and the ra- 
tio /?o - they depend practically on a single combination 
mfc-0.63m c +0.3/U 2 . The effect of the spin-orbital average 
p^ LS is negligible. The uniform dependence can be under- 
stood noting that the spectrum is mainly determined by the 
parton expression, with nonperturbative effects playing a 
relatively insignificant role. 

In spite of this 'inefficiency' of higher moments in con- 
straining heavy quark parameters, lepton moments already 
allow a decent determination of \V c bY Indeed, its value ex- 
tracted from F s i(B) has the following dependence: 

^2 = l-0.65[(m i -4.6GeV)-0.61(jn c -1.15GeV)] 
+ 0.013 (pi- 0.4 GeV 2 ) + 0.1(^,-0. 12 GeV 3 ) 
+ 0.06(/4-0.35 GeV 2 ) - 0.010>L +0.15 GeV 3 ) = 
1 - 2|| l(E e )- 1.38 GeV] - 0.06 (m c - 1.15 GeV) 

- 0.07(/4-0.4 GeV 2 ) - 0.050^-0. 12 GeV 3 ) - 
0.0804-0.35 GeV 2 ) - 0.005 (p{ s +0.15 GeV 3 ); (2) 

a combination of the parameters has been replaced by the 
first lepton moment in Eq. Q, and the residual sensitivity 
to p 2 G and p 3 LS is illustrated. The combination determining 
the semileptonic width is also very close! We see that the 
precise value of charm mass is irrelevant, but reasonable 
accuracy in p 2 and p 3 D is required. 

The first hadronic moment takes the form 

<M 2 > = 4.54 GeV 2 - 5.0 [(m/,-4.6 GeV) - 

0.62 (m c - 1.15 GeV)] 
- 0.66 Ou 2 -0.4GeV 2 ) + (p^-0.12 GeV 3 ) , (3) 
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i.e., again given by nearly the same combination m b - 
0.7m c + 0.1/4 - Q.2p 3 D as the lepton moment. Not very 
constraining, this provides, however a highly nontriv- 
ial check of the HQ expansion. For example, taking 
the DELPHI'S central value for (M 2 ,) we would predict 
(E() = 1.377 GeV, while experimentally they obtain 
{E r ) = (1.383 + 0.015) GeV [[H|. In this respect such 
a comparison is more critical than among the lepton mo- 
ments themselves. In particular, these two first moments 
together verify the heavy quark sum rule for Mj-mj, with 
the accuracy about 40 MeV! 

The dependence on heavy quark parameters expectedly 
changes for higher hadronic moments: 

{{Ml-{M 2 x )f) = 1.2 GeV 4 -0.003(m fe -4.6 GeV) 

- 0.68 (m c - 1 . 15 GeV) + 4.5 (^-0.4 GeV 2 ) 
-5.5(p^-0.12GeV 3 ) , 
<(M 2 -<M 2 » 3 > = 4 GeV 6 + (mi-4.6GeV) 

- 3 (m c -1.15 GeV) + 5 Gu 2 -0.4GeV 2 ) 

+ 13(p^-0.12GeV 3 ) . (4) 

Ideally, they would measure the kinetic and Darwin expec- 
tation values separately. At present, however, we have only 
an approximate evaluation and informative upper bound on 
p 3 D . The current sensitivity to ji 2 and fP D is about 0.1 GeV 2 
and 0.1 GeV 3 , respectively. 

The experimental constraint on the combination driving 
T s i(B) in the reported DELPHI measurements appears 
stronger for the hadronic moment. Using it instead of (Ef) 
we would arrive at 

1 + 0.14[<M 2 >-4.54GeV 2 ] - 0.03 (m c - 1.15 GeV) 
+ 0.1(ju;-0.4GeV 2 ) + 0.1(pJ,-0.12GeV 3 ). (5) 

We see that measuring the second and third hadronic mo- 
ments is an essential step in implementing the comprehen- 
sive program of extracting \V c b\ (see figures in M. Calvi's 
talk, [ 143, and in r 1161 ). Neglecting possible theoretical 
uncertainties in the above relations, we get, for example, 

\V ch \ = 0.0421 (l ± 0.01 S Lwidth + 0.015 HQp ar) (6) 

from only DELPHI hadronic moments. Incorporating into 
the fit the full set of moments they arrived at even smaller 
error interval. 

Clearly, more work - both theoretical and experimental - is 
required to fully use the potential of inclusive semileptonic 
decays. It is crucial that this extraction carries no hidden 
assumptions, and at no point we rely on 1 jm c expansion. 
Charm quark could be either heavy, or light as strange or 
up quark, without deteriorating - and rather improving - 
the accuracy. 



A similar analysis can be applied to the moments with a 
cut on lepton energy; in particular, CLEO has measured 
a few lepton energy moments for Eg > L5 GeV with un- 
precedented accuracy. As mentioned above, they also fix 
more or less the same combination of masses and nonper- 
turbative parameters. The value comes out close, but does 
not literally coincide with that obtained by DELPHI. This 
clearly deserves further scrutiny. However, since the ac- 
curacy of theoretical expressions is limited, in particular 
for relatively high cut employed, at the moment I do not 
see reasons to worry. As we heard at this Workshop N12I . 
CLEO was able to do similar measurements lowering the 
cut down to about 1 GeV. We are looking forward to the 
new data - they can be treated more reliably by theory and 
are expected to provide critical checks. 

Due to space and time limitations I have to omit discussion 
of an important question of convergence of power expan- 
sion we employ for the analysis of inclusive decays, and 
closely related to it theoretically problem of possible vio- 
lations of quark-hadron duality. Referring the reader to the 
dedicated publications [8 6 1 (an interesting complemen- 
tary discussion can also be found in the recent CKM pro- 
ceedings [ 17 1, Chapter 3, Sect. 2.3), I only mention that 
the answer radically depends on what concretely is stud- 
ied, and how. Using the proper approach the effects are 
under good control for total semileptonic widths, but not 
necessary so for higher hadronic moments where some im- 
provements may be required to achieve desirable precision. 

Experimental cuts and hardness. There is a problem, 
however, which should not be underestimated. The intrin- 
sic 'hardness' of the moments deteriorates when the cut 
on E( is imposed. As a result, say the extraordinary exper- 
imental accuracy of CLEO's R0-R2 cannot be even nearly 
utilized by theory, whether or not the expressions we use 
make this explicit. 2 

For total widths the effective energy scale parameter is gen- 
erally Q = mi,-m c . Where OPE applies we can go beyond 
purely qualitative speculations about hardness. Then it is 
typically given by Q < <w max , with <w max the threshold en- 
ergy at which the decay process kinematically disappears 
once mi, is replaced by nib-u). With the E[ >E m i n cut then 

Q^m b - E min - y]E 2 min + m 2 (7) 

constituting only meager 1 .25 GeV for E m [ n = 1.5 GeV, and 
falls even below 1 GeV for the decays with Eg > 1 .7 GeV. 
A closer look reveals that such a limitation appears relevant 
for understanding the 'unexpected' behavior of the first 
hadronic moment with respect to the cut on Eg reported 
by BaBar last summer I T19I . 

2 An instructive example of how naive analysis can miss such effects was 
given in 1 1181 . Sect. 5. 
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In b — > s + y decays one has Q m/, - 2£ m j n , once again 
a rather soft scale 1.2 GeV if the lower cut is placed at 
E y = 2 GeV. Hence, the reliability of theory can be ques- 
tioned when one aims at maximum precision. For higher 
moments the hardness further deteriorates in either decays. 
A high premium should then be put on lowering the cuts [ 

ED. 

On the theoretical side, the higher hadronic moments can 
be affected by nonperturbative physics formally scaling as 
powers of I /nib greater than 3. At the same time, these mo- 
ments are instrumental for truly model-independent com- 
prehensive studies of B mesons; improvement is needed 
already for the third moment, its expression given above is 
not too accurate. Considering alternative kinematic vari- 
ables will help to improve the convergence. This is related 
to the peculiarity of the kinematic definition of the invari- 
ant mass M x : 

M x = (P B -qf=p 2 c + 2(M B -m b )(m h -qo) + (M B -m b ) 2 (8) 

where p c = nib v—q is the c-quark momentum as it emerges 
at the quark level and nib-qo = E c is its energy; as always 
q is the lepton pair four-momentum. In the standard OPE 
we compute separately A=p 2 -m 2 having the meaning of 
the final state quark virtuality, and E c . Power expansion 
in terms of local heavy quark operators applies to both of 
them, or any their product. Convergence of A or its powers 
is satisfactory; large corrections in higher moments emerge 
from the product 2{M B -nib)E c . Moreover, it is powers of 
A that are directly related to the higher-dimension expecta- 
tion values. They also determine the distribution over the 
off-shellness of the final-state c quark, much in the same 
way as the moments of the light-cone heavy quark distribu- 
tion function shape the photon energy spectrum in b — > s+y. 

Therefore, it is advantageous to trade the traditional 
hadronic mass M x for the observable more closely corre- 
sponding to the quark virtuality, defined as 

N 2 =M 2 X -2AE X , (9) 

where Ex = M B - qo is the total hadronic energy in the 
B restframe, and A a fixed mass parameter (a constant in 
N\ does not affect higher moments). Preferred values are 
about A a Ms-nibi 1 GeV) and can be taken 500-700 MeV. 
The higher moments ((N 2 - (N 2 )) 2 ), {{N 2 X - (N 2 )) 3 }... 
should enjoy better theoretical stability. It may be possi- 
ble to further improve it dividing N x by certain kinematic 
factors. 

The kinematic variable N\ is not well constrained inclu- 
sively at LEP experiments. Is it possible to measure it at 
B-factories? This question was raised at ICHEP2002, and 
it seems the answer is positive; I'm grateful to colleagues 



from BaBar for clarifying this during the V x b and V tx Work- 
shop at SLAC. This possibility should be explored. 

Problems for HQE? So far I was seemingly quite opti- 
mistic confronting OPE-based heavy quark expansion for 
inclusive semileptonic decays to experiment. There are two 
problems, however which are often mentioned as constitut- 
ing problems for the OPE. One is a stronger than predicted 
r i20l dependence of (M x ) on the lepton cut Ef reported last 
summer by BaBar N13I . The second is a possible deviation 
of the heavy quark parameters from b — > s + y spectrum 
reportedly deteriorating the global fit when these data are 
included rl2"0l 

I am not sure we really witness here the problem for the 
theory based on the OPE. Clearly, future confirmation of 
experimental results would be necessary to draw any rad- 
ical conclusion, as has been emphasized in a number of 
talks. Yet, I think even the present data taken at face value 
do not indicate problems with applicability of the OPE, but 
may rather be attributed to certain unjustified simplifica- 
tions made in the theoretical analyses. 

The corresponding problems with the b — > s + y spectrum 
actually had been discussed in paper [ 18 1 contributed to 
the CKM 2002 Workshop. Sect. 5 specifically pointed out 
the complications arising when the cut on E y is not low 
enough. In this case the routinely used expressions for the 
moments in the form of simple sum of perturbative and 
nonperturbative terms, are not justified. This is another 
side of the same deteriorating 'hardness' I discussed ear- 
lier. Theory would become more reliable if the cut could 
be lowered by even 100 MeV. 

The similar problem plagues usefulness of the 'global fits' 
like one reported in Ref. [ 1201 - they combine accurate 
and well justified theoretical expressions with those of an 
imprecise nature and treat them on parallel footing. The 
literal outcome like the overall quality of the fit may then 
be not too meaningful, and some conclusions may happen 
to be misleading. The employed approach to treat theoret- 
ical uncertainty was also criticized by some experimental 
colleagues. While such fits are a reasonable starting point 
when new data are just coming, clearly a more thoughtful 
and sophisticated way to analyze the whole set of data has 
to be elaborated to confront data and theory seriously and 
to target the ultimate precision. 

Here I would like to dwell in more detail on the first alleged 
problem, the £>-cut dependence of (M x ), since it attracted 
much attention and caused various speculations. There 
may be a certain pure experimental controversy which will 
eventually lead to some change in numbers. Its concise dis- 
cussion has been given by M. Luke r 1211 . I abstract here 
from such possibilities and accept the data at face value. 

On one hand, I have mentioned that, say comparing the 
DELPHI'S leptonic and hadronic moments one gets an 
impressive agreement at the nonperturbative level. On 



Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 




0.8 1.0 1.2 1.4 1.6 



Figure 1. (Mf ) vs. lower cut on E ( as measured by BaBar [ 
19 1. Nearly straight line shows their quoted theoretical prediction 
given in Ref. 1 1201 . 

the other hand, the experimental measurement of (M|) at 
different E{ performed by BaBar apparently shows a far 
stronger dependence than the curve quoted as a theoretical 
prediction, Fig. 1. This may look as a controversy: do we 
observe a triumph or failure of the OPE? 

I would advocate the opinion that the discrepancy - which 
may look serious - in fact is hardly established since the 
theoretical uncertainties in the evaluation of (M^) have 
been essentially underestimated here, in particular at the 
higher end of £>. (Claimed theory error bars ranged from 
+0.035 GeV 2 to ±0.05 GeV 2 for low and high E t , respec- 
tively rl20l/) 

To give a sense of numbers without submerging into tech- 
nicalities, I recall a few relevant facts. 

• OPE computes in first instance not the plotted {M\ - 
Mjj) =i 0.4 GeV 2 for which the relative discrepancy looks 

of order unity, but rather <M 2 -(m c + A) 2 ) =s 1 .2 GeV 2 - and 
all of it is the effect of strong interactions, mostly nonper- 
turbative. 

• The actual sensitivity of (M^) to heavy quark parame- 
ters is illustrated by the fact that variation of nib by only 
±20MeV changes <M 2 > by +0.1 GeV 2 , the scale of the al- 
leged discrepancy. 3 

Rephrasing Ben Bradley, one may then suspect that the 
problem originates from "Dealing in the expressions, not 
necessarily in (OPE) truths". 

The sensitivity to the precise value of mj refers to varying 
only Mb while keeping all other heavy quark parameters 
we use, in particular m c and yu 2 fixed. Since, as mentioned 
earlier, charm mass expansion has a questionable accuracy 
and this is the only independent way to get a direct con- 

3 Z. Ligeti claimed he and collaborators did not find such a strong depen- 
dence in the theoretical expressions implying that this should undermine 
the proposed explanation. I insist that this dependence is true. 
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straint on nib - m c , this is an appropriate yardstick; even 
having m c known exactly we a priori cannot pinpoint the 
value of nib with better than 30 to 40 Me V. 

Of course, this example helps only to illustrate the scale 
of the involved theoretical uncertainties in more familiar 
terms. By itself varying nib by such a small amount can 
shift overall theoretical prediction up or down, but cannot 
noticeably change the variation of {MV) with increase of 
the E( cut up to 1 .5 GeV. Whatever nib is, it is a constant, 
the actual b quark mass. This might suggest the above con- 
sideration is not relevant for reconciling theory with the 
data. Such a conclusion, however would tacitly assume 
that other theoretical uncertainties are likewise indepen- 
dent on the placement of the cut. This would be grossly 
wrong. In fact, we know that the OPE expansion deterio- 
rates fast when E( increases and must completely blow up 
when Ei approaches the borderline at a fraction of GeV 
below 2.15 GeV. 

The numbers following Eq. @ characterizing the hard- 
ness qualitatively suggest the breakdown scale can be just 
around 1 .7 GeV. It can be concluded in a more quantita- 
tive fashion that the usual power expansion of (M 2 ) loses 
sense at least for cuts above just Ef = 1.7 GeV. To il- 
lustrate this, let us assume the chromomagnetic operator 
vanishes, and neglect the perturbative effects. The lead- 
ing term in the OPE is then given by the kinetic average 
ju 2 , and it contributes with the negative sign to (M^-m 2 ). 
There is, however a positive leading parton-like contribu- 
tion 2(MB-mb)(E c ). With the canonic values of parameters, 
mt, — 4.6 GeV and /i 2 = 0.4 GeV 2 we would obtain reason- 
able positive (M^)-M? at low lepton cut, it would vanish 
for the cut near 1.7 GeV and, retaining only the leading 
terms (M^) becomes smaller than for higher cuts (the 
shift in the meson masses when playing with // 2 and 
should be accounted for). Clearly this would be physically 
meaningless, signifying that higher order corrections must 
become 100% important, whatever are our naive estimates 
of the next-to-leading terms. 

This qualitative behavior of the OPE is quite natural. The 
question then arises - why the increasing uncertainties 
were not revealed in the corresponding theoretical calcu- 
lations? This is related to a certain subtlety in applying the 
OPE which was missed already in the original estimates [ 
1221 . Those early evaluations, representing a natural first 
step in the analysis, are conceptually incomplete in this re- 
spect. The underlying reason was addressed in Ref. Ml 81 . 
and I briefly reiterate it here. 

Inb^ctv decays this is somewhat obscured by three -body 
kinematics with massive c quark. It is more transparent for 
b — > s + y which we then can look at first. Let us consider a 
constrained fraction of the B—>X s +y events 

M B 

i r— dr fcjy 

l-O r (£) = — dEy-HT- (10) 

l bsy JE at ^y 
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Ignoring the perturbative corrections, the usual 1/m^ ex- 
pansion always yields the spectrum in the form of expand- 
ing around the free-quark kinematics: 



1 dT, 



bsy 



= a 5{E y - + Z> 5'(E y -"-±) + c 8"{E y - f ) + . 



where a, b, ... are given by the expectation values of local 
b-quark operators over B. Naively computing 1 - <f> y (E), 
or spectral moments over the restricted domain in this way 
would yield unity for any E > ^ - the result clearly un- 
justified. The actual behavior of the rates is described by 
the heavy quark distribution function. Its tail is indeed 
exponentially suppressed by a typical factor e~ cQ ^ t " mS ' at 
Q(E y ) » jUhadr, however for Q(E y ) ~ yt/hadr the true integral 
differs by terms of order 1 . 

This is an expected behavior. The point is rather that this 
is totally missed in the naive OPE and in the way to gauge 
the theoretical uncertainty based on it! Conceptually this is 
related, as illustrated below, to the limited range of conver- 
gence of the OPE for the width, determined in this case by 
the support of the heavy quark distribution function. 

The OPE for the semileptonic decays is similar in this re- 
spect. Both the lepton spectrum itself and the moments of, 
say M 2 , have similar end-point singularities in the power 
expansion, their resummation likewise yields a certain dis- 
tribution function strongly affecting moments at Q < fJ-hmb- 
In b— > u Iv this is the domain 2E( > m/,— /ihadr- 

The situation turns out more deceptive for b — > civ. In 
contrast to b — > u, charm mass leads to an apparent soften- 
ing of the end-point behavior where some of 6(Ee-E m!ix ), 
6'{Ef-E m . dx ), ... are replaced by regular functions. The lat- 
ter, in fact represent the finite-width realization of the same 
distributions, with the effective width ~ m 2 /2mb- If m 2 
were much larger than /ihadi-Wi, this would indeed damp the 
end-point nonperturbative corrections. In practice, how- 
ever, it appears that qualitatively m 2 < 2/ihadr'«fe and the 
effect of the nonperturbative smearing may still dominate. 
Then it is not reflected in the size of the naive OPE terms 
for the decay distributions. 

On the other hand, here in b — > civ decays these singu- 
lar Wilson coefficients do depend on the lepton cut through 
the 'tails' of the end-point spikes. In contrast to b— > s + y 
where total absence of the E^t-dependence alarms one that 
something is missing, the residual dependence in theoret- 
ical expressions for terms on E cut masks this, and 
their suggested size strongly underestimates the actual cor- 
rections in the upper part of the spectrum. 

This becomes rather obvious in the limit where m 2 <sc 
2jUhadrffJi [0. Here the strong nonperturbative effects com- 
pletely dominate the rates or moments at E{ > ^ m " ~ ^T"> 
while the naive OPE terms seem to suggest the correspond- 
ing effects are still suppressed by powers of m c jm\, there 



since are given by the tails of the 5-like distributions local- 
ized in the domain of m 2 /2ntt near the end point. 

This understanding reveals the proper approach to judging 
the convergence, as well as the actual uncertainty of the 
OPE series for the decay distributions. This follows from 
the consistent way of deriving the OPE for them [ 3 1: one 
should analyze the power expansion for the whole forward 
transition amplitude in question, rather than the expansion 
obtained for its absorptive part one routinely is interested 
in. 

Say, for fc->s + ywe would need to take T{qo; 0), expand 
it in 1 / {nib — 2qo) and study convergence of the power ex- 
pansion for all |m/,-2gol > nib—2E y mx . For the semileptonic 
width we would select [ 8 1 



,}</ [q 2 h(q ; q 2 ) + ^{2E t {q -E e )-^) h 2 (q {) ; q 2 ) 

+ q 2 (2E f -q )h 3 (qo-q 2 )] , (11) 



etc. Computing the width directly amounts to integrating 
these hadronic functions over a closed contour, which may 
lead to certain kinematic cancellations of individual coeffi- 
cients even if the overall result is not convergent at all. A 
simple illustration is provided by the toy function 



7X<?o;0) 



^i-^hadr-2^0 



(12) 



for which 



- J^qoJmT(q -,0)=&(^f^-E min ), (13) 

while its naive power expansion in 1 jm,b yields 
i JT d<7o Im T(q ; 0) = J dq Q 6(q -'%) 

+ dq d'(qo-f) + ... = $(f-E mm ), (14) 



i.e. exactly 1 for any cut below w/,/2. (The actual b — > 
s + y is not much different - one just needs to integrate 
this over yUh a d r with the light-cone distribution function 
F(//hadr)-) The flaw in relying on expansion dl4> is obvi- 
ous: this is not an expansion in yUhadr/'«fc, but in powers of 
jUhadr /(mb-2qu) (cf. with energy release), and it can be used 
only at 1 / (m\,-2q§) below the radius of convergence of the 
function T(qo; 0) itself, which is just mi,-2qo > jUhadr- The 
power series for Jdqolm T(qo; 0) which consists of the sin- 
gle, first term, of course converge for any qo, but may yield 
the wrong answer. 



4 Convergence for hadronic r decay widths also should be analyzed in the 
similar way. 
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Let me parenthetically note that the degree of convergence 
with a cut in E( can be also assessed by studying the OPE 
for moments containing additional powers of E[-E max , say 
{{E-E mm ) k M 2 x ) beyond just <M 2 ) we are directly interested 
in (here E miX is the end-point energy). However, this is a 
less straightforward and somewhat more ambiguous way. 

Returning to our problem, we conclude that the intrinsic 
uncertainties of the straightforward theoretical expressions 
were clearly underestimated in [ 20 1, in particular for sig- 
nificant cuts in E(. 

It is more difficult to quantify this. We know that (Mj^-M?) 
must end up at a low value -0.3 GeV 2 when one ap- 

M 2 - M 2 

proaches the endpoint Emax = £ M D smce there or >ly D 
can be produced, and (M x ) = M 2 D . Since actual {M x ) 
is a smooth function of E cnt , we would a priori expect a 
behavior similar to the curve given by BaBar, rather than 
what they showed as the theoretical prediction, at least in 
the right-most part of the plot. However, at present it is 
difficult to evaluate the actual difference: the steep fall- 
off may start around 1.7 GeV or even a little higher. Or 
a more gradual decrease over a wider domain can show 
up for E{ even below 1 .5 GeV, depending on the actual 
behavior of higher-order power corrections. Likewise it 
is hard to quantify, say the numerical difference between 
(M x )E ml =cr(M x )E CM =i.5GeV and its value as given by the first 
three terms in its naive l/m^ expansion. It can naturally be 
as large as 0.15 GeV 2 , but may turn out both smaller or 
even larger. 

The practical recommendation can be suggested as stem- 
ming from this analysis. To extract heavy quark param- 
eters or \V c b\ experiments should use as low cut on Ef as 
possible in their data sets where they are still sufficiently 
accurate and reliable. It is then instructive to compare the 
theoretically predicted (M x ) at higher cuts on £> based on 
these robust low-£> determinations, with direct experimen- 
tal data. Since the overall shift is very sensitive already 
to quark masses themselves, and the data may be strongly 
correlated, the particular quantity of interest is 



AM 2 (£ cut ) = (M 2 )^ -<M 2 > £tu , . 



(15) 



Theory must properly describe its slope at lower E cat , but 
possibly not approaching 1 .5 GeV. Confronting AM X with 
theoretical predictions in a wider range we may infer a 
nontrivial qualitative information about strong dynamics of 
heavy quarks, similar to studying the photon spectrum in 
b— > s + y. 

BPS limit. An intriguing theoretical environment opens 
up if yU 2 (l GeV) is eventually confirmed to be close enough 
to ju 2 (1 GeV) as currently suggested by experiment, say it 
does not exceed 0.45 GeV 2 . If p 2 -p 2 c <K p 2 it is advan- 
tageous to analyze strong dynamics expanding around the 



point p 2 = p? G [0. This is not just one point of a contin- 
uum in the parameter space, but a quite special 'BPS' limit 
where the heavy flavor ground state satisfies functional re- 
lations a-n\B) = 0. This limit is remarkable in many re- 
spects, for example, saturates the bound [ 23 1 g 2 > | for 
the slope of the IW function. I hasten to recall that already 
quite some time ago there were dedicated QCD sum rules 
estimates of both g 2 [ 24 1 and p 2 [ 25 1 yielding literal val- 
ues nearly at the respective lower bounds, supporting this 
limit. 

The SV heavy quark sum rules place a number of impor- 
tant constraints on the nonperturbative parameters. For in- 
stance, they yield a bound on the IW slope 



nl-H 2 c = 3s 2 (g 2 -\), 0.4GeV<e<lGeV 



(16) 



so that g 2 can barely reach 1 being rather closer to 0.85. It 
is interesting that this prediction [ 7 1 turned out in a good 
agreement with the recent lattice calculation [ f26l g 2 - 
0.83^ Jj^gj. This leaves only a small room for the slope 
of the actual B — > D* formfactor, excluding values of g 2 
in excess of 1.15 - 1.2. This would be a very constraining 
result for a number of experimental studies, in particular 
for extrapolating the B — > D* rate to zero recoil. Since 
there is a strong correlation between the extrapolated rate 
and the slope, this may change the extracted value of | V c aI- 
Therefore, it is advantageous to analyze the B — > D*(v data 
including the above constraint as an option, and I suggest 
experimental colleagues to explore this in future analyses. 

The experimentally measured slope p 2 differs from g 2 by 
heavy quark symmetry-violating corrections. The estimate 
by Neubert that p 2 is smaller than g 2 , p 2 g 2 -0.09 seems 
to be ruled out by experiment. It is not clear if a better 
estimate can be made in a trustworthy way. 

The whole set of the heavy quark sum rules is even more 
interesting. Their constraining power depends strongly on 
the actual value of p 2 . When it is at the lower end of the 
allowed interval, the BPS expansion appears the most ef- 
fective way to analyze them. 



Miracles of the BPS limit 



for nonperturbative parameters hold in the limit p 2 - 1,2 • 



A number of useful relation 

they include g 2 — |, A = 2E, p 3 LS = -p 3 D , relations for 
nonlocal correlators p\ G = -2f? m , p\ + p\ c = -(p^+p|), 
etc. 

This limit also extends a number of the heavy flavor sym- 
metry relations for the ground-state mesons to all orders in 
1/m: 

• There are no formal power corrections to the relation 
Mp-niQ + A and, therefore to m\,-m c - Mb- Mi,. 

• For the B—>D amplitude the HQ limit relation between 
the two formfactors 



2 M B -M D 2 

f - (q } = -m7TW d U{q ) 



(17) 
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does not receive power corrections. 

• For the zero-recoil B — > D amplitude all 6\/ m k terms 
vanish. 

• For the zero-recoil amplitude /+ with massless leptons 



U{{Mb-M d Y) = 



2 M B + M D 



2 y/M B M D 



(18) 



to all orders in 1 jniQ. 

• At arbitrary velocity power corrections in B — > D vanish, 



M B + M D 



. Ml+Mj-q 

S| 2M B M D 



(19) 



2 y/M B M D 

D decay rate directly yields Isgur-Wise 



so that the B 
function %(w). 

It is interesting that experimentally the slope of the B—>D 
amplitude is indeed smaller centering around p 2 D} » 1.15 [ 
IT71 indicating qualitative agreement with the BPS regime. 

What about the B — > D* amplitude, are the corrections 
suppressed as well? Unfortunately, the answer is negative. 
The structure of power corrections indeed simplifies in the 
BPS limit, however 5i/ m z, £i/ m 3 are still very significant [ 
\7\, and the literal estimate for Fd*(0) falls even below 0.9. 
Likewise, we expect too significant corrections to the shape 
of the B — » D* formfactors. Heavy quark spin symmetry 
controlling these transitions seems to be violently affected 
by strong interactions for charm. 

A physical clarification must be added at this point. Ab- 
sence of all power corrections in 1 / mq for certain relations 
may be naively interpreted as implying that they would 
hold for arbitrary, even small quark masses, say in B — > K 
transitions. This is not correct, though, for the statement 
refers only to a particular fixed order in 1 /mg expansion in 
the strict BPS limit. In fact the relations become more and 
more accurate approaching this limit only above a certain 
mass scale of order Aqcd, while below it their violation is 
of order unity regardless of proximity of the heavy quark 
ground state to BPS. 

Quantifying deviations from BPS. The BPS limit cannot 
be exact in actual QCD, and it is important to understand 
the scale of violations of its predictions. The dimensionless 
parameter describing the deviation from BPS is 



B = \\no\cm \B)\\ = ^3(q 2 -±) = 3 (£ It*"//) 5 . (20) 

Numerically B is not too small, similar in size to generic 
1 /m c expansion parameter, and the relations violated to or- 
der B may in practice be more of a qualitative nature only. 
However, fi 2 -^ 2 -, oc B 2 can be good enough. Moreover, we 
can count together powers of 1 /m c and B to judge the real 
quality of a particular HQ relation. Therefore understand- 
ing at which order in B the BPS relations get corrections is 



required. In fact, we need classification in powers of B to 
all orders in 1 jrriQ. 

Relations dl7> and dl9l for the B — > D amplitudes at arbi- 
trary velocity can get first order corrections in B. Thus they 
may be not very accurate. The same refers to equality of 
p\ c and —1p\ n . The other relations mentioned for heavy 
quark parameter are accurate up to order B 2 . The other im- 
portant BPS relations hold up to order B 2 as well: 

• Mb-Mq = m,b—m c and Mo - m c +A 

• Zero recoil matrix element {D\cjob\B) is unity up to 
0(B 2 ) 

• Experimentally measured B-^D formfactor /+ near zero 
recoil receives only second-order corrections in B to all or- 
ders in 1/mg: 



v ; 2 ^M B M D 



(21) 



This is an analogue of the Ademollo-Gatto theorem for the 
BPS expansion. 

The similar statement then applies to /_ as well, and the 
HQ limit prediction for /_//+ must be quite accurate near 
zero recoil. It can be experimentally checked in the decays 
B^Dtv t . 

As a practical application of the results based on the BPS 
expansion, one can calculate the B — > D decay ampli- 
tude near zero recoil to use this channel for the model- 
independent extraction of \V c b\ in future high-luminosity 
experiments. For power corrections we have 

M B +M D 



U ((M B -M D ) 2 ) 



1 + 







"-) 


V * J 


\m c 


m b j 



2 y/M B M D 



We see that this indeed is of the second order in B. More- 
over, A-2S is well constrained through fi^—fig by spin sum 
rules. Including perturbative corrections (which should be 
calculated in the proper renormalization scheme respecting 
BPS regime) we arrive at the estimate 



2 ^M B M D 
M B + M D 



/ + (0) = 1.03 ±0.025 



(23) 



It is valid through order B 2 ^- accounting for all powers of 
1 jniQ to order B l . Assuming the counting rules suggested 
above this corresponds to the precision up to 1/mi, essen- 
tially better than "gold-plated" B — > D* formfactor where 
already 1 /m^ terms are large and not well known. There- 
fore, the estimate J23t must be quite accurate. In fact, the 
major source of the uncertainty seems to be perturbative 
corrections, which can be refined in the straightforward 
way compared to decade-old calculations. 

Conclusions. Experiment has entered a new era of explor- 
ing B physics at the nonperturbative level, with qualitative 
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improvement in \V c b\. The comprehensive approach will 
allow to reach a percent level of reliable accuracy in trans- 
lating r s i(Z?) to |V C 2,|. Recent experiments have set solid 
grounds for dedicated future studies at B factories. We 
have observed a nontrivial consistency between quite dif- 
ferent measurements, and between experiment and QCD- 
based theory, at the nonperturbative level. 

There are obvious lessons to infer. Experiment must strive 
to weaken the cuts in inclusive measurements used in ex- 
tracting \V c b\- Close attention should be paid to higher mo- 
ments or their special combinations, as well as exploring 
complementary kinematic observables. 

The theory of heavy quark decays is now a mature branch 
of QCD. Recent experimental studies of inclusive decays 
yielded valuable information crucial - through the compre- 
hensive application of all the elements of the heavy quark 
expansion - for a number of exclusive decays as well. This 
signifies an important new stage of the heavy quark theory, 
since only a few years ago exclusive and inclusive decays 
were often viewed as largely separated, if not as antipodes 
theory-wise. 

Generally speaking, there is ample evidence that heavy 
quark symmetry undergoes significant nonperturbative cor- 
rections for charm hadrons. However, there appears a 
class of practically relevant relations which remain robust. 
They are limited to the ground-state pseudoscalar B and 
D mesons, but do not include spin symmetry in the charm 
sector. 

The accuracy of these new relations based on the proxim- 
ity to the "BPS limit" strongly depends on the actual size 
of the kinetic expectation value in B mesons, ^(1 GeV). 
The experiment must verify it with maximal possible ac- 
curacy and fidelity, without relying on ad hoc assumptions 
often made in the past with limited data available. This 
can be performed through the inclusive decays already in 
current experiments. If its value is confirmed not to ex- 
ceed 0.45 GeV 2 , the B — > D decays can be reliably treated 
by theory, and the estimate T+(0) - 1.03 can provide a 
good alternative element in the comprehensive program of 
model-independent extraction of |V c j|. 

There are many other important consequences of the BPS 
regime. The slope of the IW function must be close to 
unity, and actually below it. The related constraints on 
the slope p 2 of the experimentally observed combination of 
— > D* formfactors should be incorporated in the fits aim- 
ing at extrapolating the rate to zero recoil. Finally, I think 
that the semileptonic decays with t lepton and v T have an 
interesting potential for both inclusive and exclusive de- 
cays at high-luminosity machines. 
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